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ABSTRACT 

Previous studies of the orbital accelerations of the high altitude balloon 
satellites, Pageos and 1963-30D have shown the existence of perturbations that 
appear to be related to solar radiation pressure but of unknown mechanism. 
The normal method of computing the radiation perturbations assumes the ef- 
fective shape of the spacecraft to be spherical but in the present paper an in- 
vestigation is undertaken to assess the perturbations that may arise when the 
satellite has an ellipsoidal shape and the radiation scattered by the spacecraft 
is no longer symmetric about the line joining the satellite and the sun. Consid- 
eration is given to both diffuse and specular reflection. The study indicates 
that a slowly precessing rotation axis might explain the anomalous accelera- 
tions found in the earlier studies. 
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RADIATION PRESSURE EFFECTS ON THE ACCELERATION 
OF HIGH ALTITUDE BALLOON SATELLITES 

INTRODUCTION 

Recent studies by Fea (Reference 1) and Fea and Smith (Reference 2) have 
shown the existence of an unexplained perturbation of the acceleration of two 
high altitude spacecraft. Both the spacecraft are balloon satellites of large area 
to mass ratio and in References 1 and 2 it was suggested that the unexplained 
acceleration might be associated with or caused by solar radiation pressure. 

Figures 1 and 2, taken from References 1 and 2, show the predicted and 
observed accelerations of Pageos (1966-56A) and Dash 2 (1963-30D). The 
difference between the computed and observed accelerations is the unexplained 
perturbation. Inspection of Figures 1 and 2 indicates several important features 
of the perturbation. Firstly, the perturbation is per ; odic ; secondly, the per- 
turbation is only present when part of the orbit is in shadow (or the amplitude 
is considerably reduced), thirdly, the amplitude of the perturbation is compar- 
able to the perturbation by solar radiation pressure, and fourthly, the period of 
the perturbation is decreasing. 

There is evidence (References 3 and 4) that one of the satellites showing 
this anomalous acceleration (Pageos) is no longer spherical and that it is 
probably shaped like a prolate spheroid. If this is true, the major assumption 
made in calculating the radiation pressure perturbations, namely, that the 
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solar radiation scattered by the satellite is symmetrical about the satellite-sun 
line, no longer holds. In such circumstances it must be expected that additional 
perturbations of the orbit will arise. 

In the present paper the perturbations to the semi -major axis of the orbit 
of a satellite having the shape of a prolate spheroid are developed. The satel- 
lite is assumed to be rotating about the major or a minor axis of the spheroid 
with a period considerably less than the period of revolution of the satellite 
about the Earth. Sunlight scattered by the satellite is assumed to be reflected 
both diffusely (according to Lambert's law) and specularly. 

INCIDENT AND REFLECTED RADIATION 

Let the satellite have the shape of a prolate spheroid whose surface is 
described by the equation 

(x 2 + y 2 ) 
b o 2 

where a c is the semi -major axis (polar radius) 

b 0 is the semi-minor axis (equatorial radius) 
and the origin of the coordinate system is at the center of the spheroid with the 
z - axis corresponding to the polar radius and the x and y axes lying in the 
equator. Let the angle between the z - axis and the sun be &, then the shape of 
the cross-section normal to the sun- satellite line is an ellipse of area A(6) where 

A(<?) = rr b 0 d (2) 
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where 


1 _ sin 2 9 cos 2 tf 

» 2 = »o 2 * W 


( 3 ) 


If the satellite is spinning about its major axis then the average cross-sectional 
area over one revolution of the satellite is A(0)and the incident solar flux, .F , 
is given by 

F, = A(0)| (4) 

where S is the solar constant in erg cm” 2 sec" 1 
c is the velocity of light in cm sec" 1 

t 

If, however, the satellite is rotating about a minor axis making an angle 9 * with 
the sun-satellite line then A(£) i^ a function of time and we need its average 
value. 

Let A( £ ' ) be the average value of A (9) then 


A (9‘) 


_1_ 
2 n 



A (9) dco 


(5) 


where ^ is the angle between the sun-rotation axis plane and the major axis of 
the satellite-rotation axis plane (see Figure 3). From Figure 3 we obtain 

cos 9 = sin 9' cos co (6) 
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and hence 


A (&') = - 3 - [d] da> 


r 


(7) 


Substituting for d from Equation 3 and for 0 from Equation 6 leads to 


a A b„ 


f 

“0 "0 I 

**'> ; -J, F 


dto 


k cos 2 a)) 1 /2 


( 8 ) 


where 


■ (-S) 


sin 2 0' 


(9) 


The solution to Equation ( 8 ) is a hypergeometric function, and can be written as 


A(0') = w a 0 b 0 F^-~. 7 . 1. k ) 


(10) 


where 


n « 0 


(ii) 
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The average incident solar flux on a prolate spheroid rotating about a minor 
axis can therefore be written 

F, - A(0')f (12) 

If the satellite were spherical the solar radiation that is reflected specularly 
would be distributed evenly over the entire unit sphere surrounding the satellite. 
If, however, the satellite is prolate or oblate there will be a direction about which 
the specular reflection will be largely symmetric and which will be the effective 
direction of any forces arising from the specular reflection. This direction will, 
for reasons of symmetry, lie in the plane containing the rotation axis of the 
satellite and the sun. We now make the first major assumption; that the effective 
direction of reflection of the specular flux is determined by Snell's law on the 
incident ray that passes through the center of the satellite (see Figure 4). We 
also make the assumption that the magnitude of the flux reflected in this direction 
approximates to that which would be reflected by a sphere of surface area equal 
to that of the spheroid. Hence the specularly reflected flux, E s , can be written 

E s = a , •Jr(t) < 13 ) 

where 

A = A(£) for rotation about the major axis 

A - A(0 ' ) for rotation about a minor axis 
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and a s is the specular albedo of the satellite. 

The assumption concerning the direction of the reflection holds for p = 0, 

V2 and ti (also f > ' ), and for 0 < 3 ^ / 2 and tt/ 2 <3 : n the direction of reflection is 
moved towards the minor axis direction as indicated by Snell’s law. Thus the 
assumption is considered adequate for the present study. 

Similar arguments can be applied to the magnitude of the reflection; the 
magnitude of the flux for 3 = 0 is probably overestimated but underestimated 
for 3 = V2. An exact formulation of the magnitude and effective direction of 
the specular reflection is extremely complex and is, at present, believed to be 
unnecessary for the present study. 

Diffusely reflected radiation is normally symmetric about the no. lui to the 
surface and this is the assumption made here (see Figure 4). The satellite is 
assumed to be a uniformly diffuse reflector (Lambert’s law) and for the purposes 
o' calculating the dependence of the magnitude on the phase (not the size) the 
satellite is assumed to be spherical. 

In Reference 5 the author has derived an expression for diffusely reflected 
radiation falling on a unit area distant r from the satellite. By putting r = ’ 
this expression may be used to give the flux (E D ) reflected diffusely in the direc- 
tion of the normal to the surface on the incident ray that passes through the 
center of the satellite. We therefore have 

= T a D ^2 (f) [ (77 ' cos€ + sin f ] (15) 
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where a D is the diffuse albedo of the satellite, 


* = 7 - (0 ♦ 4) 


and 


( 16 ) 


tan <t> 



for rotation about the major axis, or where 


and 


6 = e' - 4 >' 


tan 4 >* 



(17) 


for rotation about a minor axis. 

For convenience, we summarize the directions of the incident and emitted 
fluxes as follows: 

(a) The incident flux, Fj , is directed radially from the sun through the 
center of the satellite; its magnitude is given by Equations 4 or 12 ; 

(b) The specularly reflected flux, E s , lies in the plane containing the sun 
and axis of rotation of the spacecraft and makes an angle 2 e with the 
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incident ray through the center of the satellite; the magnitude is given 
by Equation 13, 

(c) The diffusely reflected flux, E D , lies in the plane containing the sun and 
the axis of rotation of the spacecraft and m xes an angle e with the 
incident ray through the center of the satellite, the magnitude is given by 
Equation 15. 

As a check on the effect of the approximations we have made we can integrate 
the specular and diffuse fluxes (Equations 13 and 15) over the entire unit sphere. 
This integration leads to 


total reflected flux 



(18) 


If (a s + a D ) = 1, then the total reflected flux is equal to the total incident fiux 
(Equation 4) which means that the approximations (on average) onl\ al ect the 
relative magnitudes of the specular and diffuse components and not the total 
flux. If (a s + a D ) < 1, the satellite is absorbing some of the incident radiation 
and implies we are making the additional assumption that when the satellite re- 
emits the absorbed radiation it does so isotropically so that there is no change 
in momentium of the satellite. 


PERTURBATIONS OF THE SEMI-MAJOR AXIS 

The semi-major aids is a measure of the energy of the orbit and hence the 

i 

perturbations to the semi-major axis ar' -^qual to tne work done by the forces 
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of radiation pressure on the satellite. The perturbations to the semi-major 
axis can therefore be written as (Reference 6) 


2F f "} shadow entry 

- — r cos <£ 

n am [ J shadow exit 


( 19 ) 


where 

A a is the change in the semi -major axis per revolution 
F is the flux of radiation (F is neg: tive) 
n is the mean motion 
a is the semi-major axis 
m is the mass of the satellite 
r is the geocentric radial distance of the satellite 
<P is the angle between the sun and the satellite (see Figure 5) 
Applying equation 19 to the incident and reflected fluxes already derived we 
obtain 


Aa 


n 2 a m 


;{ F i 


cos < p Q + E d cos <£j + E s cos <p. 


■>] 


shadow entry 


shadow exit 


( 20 ) 


whereto , and £ 2 are the angles between the satellite and the incident, diffuse 
and specularly reflected flux directions (see Figure 5). 

Subsequently, we shall want to allow the spin axis to precess about another 
axis so we shall assume we know the direction of the precession axis (w, 8 ) 
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with respect to the sun (see Figure 5), the position ol the sun (a, 5), the position 
of the spin axis (a s ,» s ) with respect to the precesson axis (see Figure 5) and 
the position of the satellite (a 0 , 8 0 ). 

With the aid of Figure 5 we can derive the right ascension (a) and declination 
( S) of the precession axis from 


sin S = cos r sin 5 + sin 9 cos b cos j 


- s in:' sin a. 

sin (a - a) = * 


( 21 ) 


and the spin axis (a s , S^) from 


sin 8 - cos 9 sin e + sin 6 cos 8 cos a- 

s s s s 

( 22 ) 

sin 9 ^ sin a; 

sin (a % - a) = 1 

S 


We can also derive the position of the spin axis (^ 0 , &) with respect to the sun 
from 


cos 9 


sin b sin 8 


sin oj 0 - 


+ cos 8 cos 5 


cos 8 t sin (a s 
sin 9 


cos (a s - 
~ a) 


a ) 


(23) 
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and the right ascensions and declinations of the diffuse (a , S t ) and specular 
( a 2 > ^j) reflections from 


sin Sj - sin S cos (6 + t) + cos 5 sin (6 + e) cos a> 0 


sin (6 + €) sinct 0 
sin (a t ~ a) = cos S l 


(24) 


sin S 2 


sin 5 cos (0 + 2e ) + cos S sin (# + 2e) cos a> 0 


sin (a.j - a) 


sin (0 + 2f) 
cos § 2 


(25) 


The above equations have been derived for rotation about the major axis of the 
sate'lite. For rotation about a minor axis we replace c? with 9' and t with -e in 
Equations 24 and 25. 

We are now in a position to determine <t> l and <t> 2 from 


cos® 1 = sin8 0 sin8j + cos S 0 cos Sj cos (a 0 - > (26) 


cos <j> 2 - sin S 0 sin 5 2 + cos £ 0 cos S 2 cos (a 0 - a 2 ) (27) 


which, together with. 


ccs <P 0 - sin sin b + cos $ 0 cos S cos (a 0 - a) 


(28) 


1] 



enables equation 20 to be evaluated. Equations 20 through 28 therefore represent 
the perturbation to the semi-major axis due to incident, diffuse and specular 
radiation. 

EXAMPLE 

The magnitude and form of the perturbation can best be demonstrated by an 
example. Let us assume we have a circular orbit with inclination 90 degrees, 
and the sun on the equator. Because the shadow entry and exit points are 
symmetric about the Earth-sun line Equation 20 reduces to 


Aa 


n 2 m 


[ E d 005 


+ E 

S 


COS 



shadow entry 
shadow exit 


(29) 


For this particular example it is preferable to use slightly different formula- 
tions for cos and cos <t> 2 » Let 4>’ 0 be the value of <p Q at the point of entry and 
exit into the shadow, then we can write (see Figure 6) 

cos <p lN - cos (<?+€) cos <Pq + sin {6 + e) sin<£ 0 ' cos (o>j + a> 2 ) (30) 

cos<£ lx = cos (0 + e) cos 0 O ' + sin (#+«) sin <£ 0 ' cos (31) 

and 


cos <*> 2N - cos(£+2e) cos<£ 0 ' + sin {9 + 2e ) sin 4> 0 ' cos (o^ + Wj) (32) 
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(33) 


cos £ 2X - cos (6 + 2e) eos*P 0 ' + sin (6 + 2e) sin$ 0 ' cos 

where <j> 1N , 4> 2H , <P 2H are the shadow entry and exit values of 4> l and 4> 2 . 

The angles cc l and ai 2 are defined in Figure 6. 

We can now write 

Aa = sin^ 0 ' sin (u^ + -y^sin ~j~ JEjj sin (0 + e) + E s sin {& + 2e)J (34) 

The sin <t> 0 ' term changes slowly as the orbit moves with respect to the sun but 
its sign is always positive. The sin u> 2 / 2 term is zero when there is no shadow 
on the orbit and the perturbation vanishes. When non-zero, the term is always 
positive. The term containing F- D and E s changes only slowly as the spin axis 
moves and is normally of constant sign. However, for an oblate satellite the 
sine terms can change sign if 2 e> d The sin (^, + <^ 2 / 2 ) is also positive 
unless ^ j is negative, implying that the spin axis lies near the equatorial plane 
and between the sun and the shadow points on the orbit. 

Let us further simplify our example by having the precession axis on the 
equator, so that co = n/2 on I = 0, and let the spin axis rotate slowly about the 
precession axis. The sin (u l + u 2 /2) term in equation 34 will then oscillate 
about zero with a period equal to the precession period and produce a quasi- 
periodic perturbation about the mean value. This situation could have existed 
for the Pageos satellite in July 1968. Reference 4 indicates that the spin axis of 
Pageos was near the equator on July 4th and if the axis was precessing about 
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a point even nearer the equator we should expect to observe a quasi -periodic 
perturbation in the semi-major axis. The predicted maximum magnitude of 
this perturbation of Pageos, assuming a specular albedo of 0.86 (Reference 4) 
and negligible diffuse albedo, a mass of 55 kg, a mean motion of 8 revolutions/ 
day is about 1 km/day or an acceleration of 8 x 10“ 4 revolutions/day 2 in mean 
anomaly. The observed amplitude of the acceleration shown in Figure 1 is about 
l x 10“ 4 revolutions/day 2 . When account is taken of the trigonometric terms in 
equation 34 and of smoothing in the observational data; the amplitudes of the 
computed and observed perturbations are about equal. 

CONCLUSIONS 

Expressions have been developed for the perturbation of the semi-major 
axis of the orbit of a satellite with elliptical cross-section due to 3olar radia- 
tion pressure when both specular and diffuse reflections are taken into account. 
The theory has been applied to a very simplified example resembling the orbit, 
of the Pageos satellite in July 1968 and it has been shown that if the spin axis 
is permitted to rotate about a fixed precession axis a perturbation of the semi- 
major axis is predicted which has approximately the same form and amplitude 
that is actually observed in the Pageos orbit. 

The foregoing theory and example do not necessarily explain the perturba- 
tions in the Pageos and 1963-30D orbits but do suggest that a mechanism of the 
type described here could be the explanation. A more detailed examination of 
the theory and its application to these two satellites is being undertaken. 
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Figure 2. 1963 30D, Observed and Computed Orbital Acceleration v. Modified Julian Date (Nov. 1963 to Jun. 1965) 

(after Fea and Smith, Reference 2) 
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Figure 3. Rotation of Satellite About a Minor Axis. 
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Figure 4. Incident and Reflected Radiation. 
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Figure 5, 


Relationship between Sun, Precession and Spin Axes, Diffuse and Specular 
Reflections, for Satellite Rotation about Major Axis. 
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Figure 6. Shadow Entry and Exit Geometry foi j Circular Orbit. 


21 


■shadow" 


